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Resonant tunneling through semimetallic ErAs quantum 
wells embedded in GaAs structures with AlAs barriers was 
recently found to exhibit an intriguing behavior in magnetic 
fields: a peak splitting occured only in fields perpendicular to 
the film while a second resonant channel opened for in-plane 
fields. The behavior is explained in terms of the valence bands 
states in ErAs in the vicinity of the T point, their exchange 
splitting induced by the localized Er 4f magnetic moments, 
and by a selection rule involving the total angular momentum 
component along the normal to the interface. 

PACS 73.40.Gk, 73.40.Vz 



The utilization of thin metallic layers as active com- 
ponents in semiconductor devices such as metal base 
and resonant tunneling transistors opens new avenues 
in micro-electronics. Resonant tunneling through thin 
metallic layers embedded in, a semiconductor has been 
reported by several authors £1 A promising materials sys- 
tem for the study of this phenomenon is semimetal- 
lic ErAs embedded in GaAs since the two are closely 
lattice matched.cl Because of the open 4f shell of Er- 
bium, these hetejepstructures also have interesting mag- 
netic behavior.Sia Here we provide a theoretical expla- 
nation for the unusual behavior of resonant tunneling 
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through ErAs layers in magnetic fields of different ori 
cntation which was recently reported by Brehmer et al 
and Zhang et al..H These experiments were carried out 
on resonant tunneling diode (RTD) structures consisting 
of ErAs quantum wells and AlAs barriers sandwiched be- 
tween the n+ doped GaAs substrate and GaAs capping 
layers. 

The measurements of differential conductance (dl/dV 
vs. V)B, reproduced in Fig. 0, indicate the presence of 
two different resonant tunneling channels called A and B. 
The peak A splits in a magnetic field perpendicular to the 
layers, but remains nearly unaffected (slightly shifted) by 
a magnetic field parallel to the layers. The B channel, on 
the other hand, shows no observable splitting for either 
direction of the field, is only weakly resolved as a shoulder 
in zero or perpendicular field, but is strongly enhanced 
by a field parallel to the film. The most extensive results 
and the ones we focus on here were obtained for the [113] 
orientation of the film. However, similar results were ob- 
tained for other orientations of the films. That indicates 



that the phenomena are not determined by specific intrin- 
sic crystallographic directions but rather by the relative 
orientation of field and interface. 

In the present letter, we will show that this behavior 
is explained by tunneling of electrons into the ErAs light 
and heavy hole states in the vicinity of the T-point for 
the A and B resonances respectively. In the absence of 
a magnetic field or in a perpendicular field, tunneling is 
predominantly only into the light hole bands because of 
an approximate conservation of the total angular (orbital 
and spin) momentum component along the normal of the 
film. The exchange splitting of these states induced by 
coupling to the open 4f shell leads to an anomalously 
large Zeeman splitting. On the other hand, a parallel 
field allows for tunneling into the heavy hole band leading 
to an enhancement of channel B. The splitting in this case 
will be shown to be reduced. 

To understand the tunneling behavior, a detailed un- 
derstanding of the ErAs band structure and magnetic 
properties is required. This is rfiffiovided by our recent 
electronic structure calculations™ of ErAs and related 
Er^Sci-^As alloys, the results of which were in-good 
agreement with magnetotransport measurements. 13 The 
first question to be addressed is whether the tunnel- 
ing takes place into the conduction (electron) or valence 
(hole) states. This is expected to depend on the film 
direction. For the [113] direction under consideration, 
the band structure (see Fig. ||) shows clearly that un- 
occupied electron states are more than 1.5 eV above the 
Fermi level. Hence they cannot be involved in the tun- 
neling. In fact, although the Schottky barrier height at 
the ErAs / AlAs interface is presently unknown, the Fermi 
level in the ErAs must line up with the donor states near 
the conduction band edge deep in the n+ doped GaAs 
region. Implicit in the above argument is that fcn is con- 
served in the resonant tunneling and that only electrons 
with fcii w are involved in the transport from GaAs 
because the conduction band minimum in GaAs is at T. 

Using the simple infinite barrier model, a first estimate 
can be obtained for the quantum confined states in the 
ErAs quantum well. For a film thickness of N monolay- 
ers (ML), there are n = 1 . . . N quantized values of wave 
vector = mr/Na_i_, where a± is the distance between 
two monolayers. The corresponding values of fcj_ are in- 
dicated as vertical lines in Fig. || for a layer thickness 
of 12 ML. This shows that the only the n = 1 quantized 
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hole states lie close to the Fermi level and are likely to 
be involved in the tunneling. 

To relate these quantized energies E n to the voltages 
V at which resonance occurs, we need to know the volt- 
age profile in the RTD. This is somewhat uncertain be- 
cause the top and bottom interfaces of ErAs are of differ- 
ent quality due to difficulties in the epitaxial overgrowth 
on ErAs. For a symmetric structure, there would be 
equal voltage drops across the two insulating GaAs space 
charge layers and AlAs barriers leading to a resonance 
when V/2 as E n . Direct measurements of the voltage 
drop on three terminal devices, however have shown that 
most of tke voltage drop occurs on the substrate to ErAs 
interfaced This implies that E n %V. 

Unlike our earlier calculations!, the band structure 
shown in Fig. || includes the effect of spin-orbit cou- 
pling which, as will be shown below, is essential for 
understanding the magnetic field dependence. Other- 
wise, the details of the calculation method are similar 
to those of Ref.tl The calculations are performed within 
the atomic sphere ..approximation of the linear muffin- 
tin orbital methodic and are based on local spin den- 
sity functional theory. The spin-orbit coupling terms are 
added as a perturbation to the spin-dependent Hamilto- 
nian, which is then diagonalized numerically without fur- 
ther approximations .113 The spin-polarization arises from 
the alignment of unpaired spins in the 4f level which is 
treated as an open shell core-like level without dispersion. 
In the absence of a magnetic field, the 4f localized spin 
magnetic moments are randomly oriented and the bands 
are described by a calculation without spin-polarization. 

The manifold of valence states around V can be. 
described by the Kohn-Luttinger (KL) HamiltonianjLil 
which accounts for the behavior up to quadratic terms 
in k. In the following, we use this Hamiltonian in the 
spherical approximation (i.e. 72 73 in terms of the 
conventional Luttinger parametersEd) in order to gain a 
qualitative understanding of the magnetic field effects. 
We stress, however, that this description of the nature 
of the quantum states involved in the tunneling and the 
associated conservation laws remain approximately valid 
even in the complete calculation. The KL-Hamiltonian 
we start from is: 

H= - (71 + 4 72 )k 2 + 6 72 (L • k) 2 

+ (A + A 1 k 2 + A 2 (L-k) 2 )(S-B) 

+ ^A so (L-S-i), (1) 

where L is the I — 1 orbital momentum operator de- 
scribing p-like hole states at the top of the valence band 
of ErAs, S is the spin operator and B a unit vector 
along the magnetic fi<^4- The first two terms are the 
usual Luttinger terms .LH The last term is the spin-orbit 
coupling term with the 1/2 inside the parenthesis be- 
ing a reference-level shift to the valence band maximum. 
The parameters Xi describe the splitting of the valence 
states in the presence of a magnetic field. Here the split- 



ting is primarily due to the exchange interaction of the 
holes with the ferromagnetically aligned localized mo- 
ments of 4f electrons. The alignment of the latter is welk 
described by standard Brillouin theory of paramagnets.El 
In our numerical calculations we assume saturated spin- 
polarization. This is valid because the measurements 
were taken at T = 4K and the applied magnetic field 
was 8 Tjj Direct coupling of the angular momentum to 
the magnetic field (the Zeeman effect) is neglected as it 
is much smaller than the coupling to the exchange field. 
Also neglected are the terms that lead to the formation 
of Landau levels. The above Hamiltonian is thus equiv- 
alent to the one used in our first-principles calculations 
and its parameters can be extracted from the latter. 

For reasons already explained above, we can now set 
fen = and choose the z axis along the normal to the 
interface. The quantization of k± in the infinite bar- 
rier model then fixes the k completely. In the absence 
of the spin-orbit interaction the above Hamiltonian (|l|) 
describes two sets of spin-polarized bands that are in- 
dependent of the direction of the magnetic field. The 
spin-orbit interaction now couples the spin-direction (de- 
termined by the magnetic field) to the orbital angular 
momentum direction (fixed by k± to be normal to the 
interface) . This is what leads to the non-trivial magnetic 
field orientation dependence. More precisely, without the 
magnetic field term but including spin-orbit interaction 
the Hamiltonian commutes with J z where J z = L z + S z 
is the projection of the total angular momentum onto the 
direction normal to the interface. If the magnetic field 
is along the same direction (B || z), the Hamiltonian (|l|) 
maintains this cylindrical symmetry and splits up into 
two 2x2 submatrices H^ ±1 ^ 2 — describing the cou- 
pling between the |3/2,±l/2) and |l/2,±l/2) states, 

U2A n ±B nT C n ) ^(2A n T2B nT C n ) \ 
£{2A n t 2B n t C„) \(AA n ±B n ±2C n )- A so ) ' 

(2) 

where E° = -(71 + ^i)k\ n , A n = 672^^, B n = 
Ao + Aifcj_ n , and C n = A 2 fc 2 _„. We will primarily be 
concerned with n = 1 states. Similarly the eigenvalues of 
the |3/2, ±3/2) states are 

E n ,± 3/ 2 =K + A n ± D n /2, (3) 

where D n = B n + C n . 

In the case of the magnetic field parallel to the 
interface, J z no longer commutes with the Hamilto- 
nian and there is an additional mixing of the eigen- 
states of the original Hamiltonian (without field). 
The energy spectrum of the quantum well can, how- 
ever, be found as the eigenvalues of the two 3 x 
3 submatrices h\ ± — E®, written in the basis 
^(|3/2,3/2)±|3/2,-3/2)) J ^(|3/2,l/2)±|3/2 J -l/2)), 
and -^(|l/2, 1/2) T |l/2, -1/2)) 



2 



/ A, 



V12 " V6 

. \{A n ±B n ) ^{2A nT B n ) 

\±-L(D n ) f(24TB„) \{AA n ±B n )- A so 



(4) 

The analysis of the structure of the matrices (^-||) al- 
lows one to explain semiquantitatively all the features of 
the resonant tunneling experiments in a magnetic field. 
First, we recall that the electrons tunneling form the 
n + GaAs conduction band minimum are s-like and have 
Mj = ±1/2. With the magnetic field absent or per- 
pendicular to the layers, Mj is a good quantum num- 
ber that must be (approximately) conserved. Hence tun- 
neling transitions between states with \Mj\ = 1/2 and 
\Mj\ =3/2 are forbidden. This explains why in that case 
only the A-channel (corresponding to tunneling into the 
light hole states with Mj = ±1/2) is pronounced. The 
small shoulder of the B-channel is due to the fact that 
in the real band-structure the axial symmetry is slightly 
broken. For parallel fields, the Mj states are mixed and 
the B-channel becomes available for tunneling. The ob- 
served splitting in the perpendicular magnetic field is also 
easily understood since the energies depend on the sign 
of Mj. The magnetic field dependence of this splitting 
was shown by Zhang et al.Q to follow the Brillouin theory 
of paramagnetic alignment of localized spins as expected 
in the present theory. 

Fig. ||(a,b) shows the eigenvalues of the above ma- 
trices in perpendicular and in-plane field respectively, 
with the parameters E n = -84.3fc^, A n = 59.7kj_, 
B n = — 9.75fcj_, C n = 14fc^, with kj_ in units n/aj_ 
and Aso — 0.45 and all energies in eV. These parame- 
ters were extracted from our first-principles calculations 
and appropriately adjusted for the spherical approxima- 
tion. The small spin-splitting at T, Ao ~r-J- meV was 
neglected. As discussed in our previous workB, the Fermi- 
level position was slightly shifted from its LDA value as 
it was found-, that the LDA overestimates the Fermi sur- 
face volumea because of the neglect of electron interac- 
tion self-energy effects. This figure clearly shows that 
the light-hole spin-splitting is about 5 times larger for 
the perpendicular than for the in-plane field. Also indi- 
cated in this figure are the experimental resonance posi- 
tions for varying thicknesses (i.e. different fcj_ jrl ) assum- 
ing V = E n . Overall, the agreement is quite good for 
the absolute energies of the resonances, the dispersions, 
and the light-hole heavy hole splitting. The spin-splitting 
appears to be slightly underestimated in the calculation. 
Experimentally the splitting was seen to be about 80 
meV, while the model gives 50 meV. 

The advantage of the present analytical model is that 
we can easily investigate the effect of varying parame- 
ters. Because spin-orbit interaction is fairly large, the 
coupling to the spin-orbit split-off bands can be treated 
in perturbation theory. Up to first order terms in 1/ Asoi 
the light hole spin splitting in perpendicular field is: 
A^ c (lh) « {B n - C„)/3 - AA n {2B n + C n )/9A so , which 



is of order 50 meV, while for the in-plane field, it is: 

Alc(lh) « 2B n /3 - 4A n B n /9Aso, which is of order 10 
meV. Similarly, we find that the heavy hole splitting for 
the in-plane field up to terms in 1/Aso is zero, which 
explains why no spin-splitting is observed for the B- 
channel. For the perpendicular field, the heavy hole split- 
ting is B n + C n w 40 meV but is not observed because 
in that case this channel is reduced to a weak shoulder 
on a rising background. Again up to the same order in 
I/A50, the component of the allowed J = 1/2 channel 
in the light hole eigenstate for in-plane field is given by 
{1 - [(D n /Vl2)(l + 2A n /A so )/(2A n /3)} 2 } 1 / 2 which is 
about 94%. This is consistent with the fact that the A- 
channel is not noticeably reduced in intensity when the 
in-plane field is switched on. 

A quantitative calculation of resonance intensity for 
the B channel due to the weak breaking of axial symme- 
try and its gradual switching on in the magnetic in-plane 
field will require a more detailed study of the quantum 
tunneling process. Finally, we note that the resonances in 
this system are relatively broad compared to typical reso- 
nant tunneling in semiconductor quantum wells. Among 
other broadening mechanisms, this may result from the 
loosening of k\\ conservation due to the roughness of the 
top layer interface. We thus anticipate that with further 
improvements in the material some of the finer details of 
our model may be checked experimentally in future work. 

We have also carried out fully first-principles calcula- 
tions of the bands in the [113] direction in parallel and 
perpendicular magnetic field including the spin-orbit cou- 
pling. These calculations qualitatively confirm the above 
results, but lead to a somewhat smaller light hole spin 
splitting (30 meV) in perpendicular field than the spher- 
ical model (which was already too small) and slightly 
lower masses and hence slightly lower level positions for 
a given layer thickness. A possible origin for this discrep- 
ancy is that our present calculations are based entirely 
on bulk band structures. In the quantum wells, mixing 
of Er 5d derived bands, which have larger spin-splitting, 
into the relevant hole states will occur and may lead to a 
substantial enhancement of the exchange splitting effects. 

In conclusion, the spherical Kohn-Luttinger model of 
the valence band maximum supplemented with a strongly 
enhanced spin Zeeman effect (due to coupling of the va- 
lence electrons with the localized magnetic moments of 
the 4f electrons which are aligned in the field) provides 
a good semiquantitative description of the ErAs valence 
band maximum and its associated eigenstates. The tun- 
neling behavior can be understood in terms of approxi- 
mate conservation of the total angular momentum com- 
ponent along the normal to the interface. For perpen- 
dicular or zero field this forbids the transitions into the 
heavy hole bands while the latter become allowed by mix- 
ing of states for an in-plane field. The origin of this 
approximate selection rule is that the hi ~ of the elec- 
trons originating in the GaAs conduction band is con- 
served. This in turn fixes the quantization axis for the 
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orbital angular momentum along the normal to the in- 
terface. The spin direction on the other hand is fixed by 
the magnetic field. The spin-orbit coupling is thus essen- 
tially responsible for the qualitatively different behavior 
of the tunneling in the parallel and perpendicular fields. 
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FIG. 1. Differential conductance of 12.8 ML ErAs resonant 
tunnelingrdiode in perpendicular and in-plane magnetic fields 
from Ref.EI. The A- and B-resonant channels are indicated. 
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FIG. 2. Energy band structure of ErAs along the [113] di- 
rection in zero field. The vertical lines indicate quantized k± in 
values in an infinite barrier quantum well of 12 ML. 
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FIG. 3. Band structure of ErAs near the valence band max- 
imum in the [113] direction in a saturating magnetic field 
which in (a) is perpendicular and in (b) is parallel to the 
(113) plane within the spherical Kohn-Luttinger model. Ex- 
perimental resonance positions given by V — E n (see text) 
for n = 1 and different film thicknesses (or k±^ n values) are 
indicated by circles for the A-channel and squares for the 
B-channel. The experimental spin-split A-channel resonances 
are indicated by triangles for one layer thickness in (a). The 
heavy hole (hh), light hole (lh) and split-off hole (sh) bands 
are labeled. 
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